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Abstract. 

At low temperatures, phonon scattering can become so weak that phonon transport becomes 
ballistic. We calculate the ballistic phonon conductance G for membranes using elasticity theory, 
considering the transition from three to two dimensions. We discuss the temperature and 
thickness dependence and especially concentrate on the issue of material parameters. For all 
membrane thicknesses, the best conductors have, counter-intuitively, the lowest speed of sound. 



1. Introduction 

Recent advances in nanoscience are expanding the limits of phononic thermal transport, both 
in the low conductance and high conductance side [1]. The highest thermal conductance of 
a particular material is achieved, when the phonons do not scatter at all; this is the ballistic 
transport limit. At room temperature the ballistic limit is quite hard to achieve, but at low 
temperatures (a few K and below) phonon-phonon and phonon-impurity scattering become so 
weak that ballistic transport can be observed [21 [3], and the emitted power has the typical 
Stefan-Boltzmann form [1] 

P = AaT'^, where for phonons a = T:^k%/{l^h^)T.ei/cj , (1) 

summing over the different phonon modes with speeds of sound q and radiator emissivities e^. 

In addition to bulk samples, ballistic phonon transport was also observed for thin amorphous 
SiNx membranes recently [5]. At low temperatures, the dominant thermal phonon wavelength 
becomes comparable to the membrane thickness (for a Ijim SiNx membrane this happens at 
wlOOmK) leading to an effectively 2D phonon gas. We have recently considered ballistic and 
diffusive thermal conduction in thin membranes [H [71 [8] , spanning the transition from the 3D to 
the fully 2D limit, using the correct low-frequency modes known from elasticity theory (Lamb 
modes). Here we explicitly discuss the dependence of ballistic thermal conduction on materials 
parameters. The simple question we would like to answer is: what kind of material is the best 
(worst) ballistic thermal conductor for any thickness of the membrane, and can we even define a 
material-only dependent thermal conductivity that does not depend on the sample dimensions? 

For the more typical cases where phonon transport is not ballistic, but limited by different 
scattering mechanisms in the sample, it is possible to define a material dependent local thermal 
conductivity k, which depends on the specific heat cv, average speed of sound c and phonon 



mean free path A as k = l/ScycA. For 3D crystals, cy is proportional to l/c"^, but because A 
scales as A oc c", where n = 3, 4, 5 for phonon-electron, impurity, and phonon-phonon scattering, 
respectively, k always grows with higher speed of sound [1]. This is the intuitive result that high 
speed of sound materials are better thermal conductors (like diamond). However, for the case 
of ballistic conduction this is not true anymore. The thermal conductance G (one cannot define 
a local K anymore) is calculated from the total power P that is radiated from an object as 
G = dP/dT. G actually has a dependence G oc as can be seen from the Stefan-Boltzmann 
law, Eq. (pP). The material with the lowest speed of sound is the best thermal conductor! In the 
following we will investigate what is the situation for thin membranes. 



2. Membrane eigenmodes 

In isotropic 3D bulk systems there are three independent phonon modes, two transversally and 
one longitudinally polarized, with sound velocities q and q, respectively. In the presence of 
boundaries, the bulk phonon modes couple to each other and form a new set of eigenmodes, 
which in the case of a free standing membrane are horizontal shear modes (h) and symmetric 
(s) and antisymmetric (a) Lamb modes [S]. The frequency uJh for the h modes is simply 
uJh = '^t^k'^\ + (mTr/d)'^, where fcy is the wave vector component parallel to membrane surfaces, 
d is the membrane thickness and the integer m is the branch number. The dispersion relations 
of the s and a modes cannot be given in a closed form, but have to be calculated numerically 
|10j . The lowest three branches, dominant for thin membranes at low temperatures, have the low 
frequency expressions uJh,o = ctkp cvg^ = Csk\\ and a;a,o = 2m* ^jj' '^ith the effective sound velocity 

-1 

of 



Cs = '^ctyj {cf — c^)/cf of the s mode and the effective mass m* = Ti 2ctd-^ {cf — c^)/3c 
the a mode "particle" [TO]. This lowest a mode has a quadratic dispersion instead of the usual 
phonon-like linear one, and is mostly responsible for the non-trivial behavior of thin membranes 
at low temperatures. 



3. Ballistic conduction 

Figure [1] shows a schematic of the system in consideration. A thin metal film of perimeter 
I heats the membrane directly below [7j. To simplify the discussion, we assume that the 
hot phonons from the heater have a thermal distribution and are radiated from its perimeter 
into the membrane, similar to black body radiation. We further assume that no radiation is 
backscattered. With these assumptions the total heat flow out of the detector is 
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where n(u>,T) is the Bose-Einstein distribution and a and m are the mode and branch indices 
If enough branches are used, the 3D to 2D transition can be computed from this expression. 
If the membrane is thin and temperature low (Td <^ hct/2kB), only the lowest branch of each 
mode is occupied (m=0), and we are fully in the 2D limit, in which case we get from Eq. ([2]) 
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Note that the effective mass of the lowest a mode depends on the membrane thickness and hence 
in the 2D limit P oc 1/^/d. In the 3D limit Td ^ hct/2kB, the dominant phonon wavelength is 
much smaller than d, leading to decoupling of the longitudinal and transversal modes and 



nHdhf 2 1 \ fksTY 



(4) 



As expected, P3D oc d. This means that at a fixed temperature the radiated power wiU first 
decrease with decreasing d, then reach a global minimum and will increase again, if we decrease 
d further. The minimum is approximately at the 2D-3D crossover thickness dc = hct/{2kBT). 




Figure 1. Schematic of the considered 
geometry. A thin metal film heats the part of 
the membrane directly below. The radiating 
area is A = I x d 



Figure 2. (a) Ballistic conductance G and 
(b) Ballistic conductance per unit area G/ld 
as function of membrane thickness d. The 
dash-dotted line in (a) is oc 1/d. 



4. Numerical results 

Using the 300 lowest modes, we have numerically computed the ballistic phonon conductance 
G for a range of parameters that spans the transition from 2D to 3D. Fig. [2] (a) shows the 
computed G for T=100 mK and 1=1.6 jim as a function of the membrane thickness d for four 
different materials: Pb (q = 1025 m/s, q = 2370 m/s), Al (q = 3280 m/s, q = 6580 m/s), SiN 
(q = 6200 m/s, Q = 10300 m/s), and poly-diamond (q = 10900 m/s, q = 18100 m/s). The 
values of c are calculated values for the isotropic case, based on measured single crystal elastic 
constants [12j . All materials have a physically relevant critical thickness, ranging from dc = 40 
nm for Pb to dc = 400 nm in diamond, below which G increases with decreasing thickness. The 
material with highest G is Pb, as it has the lowest speeds of sound, and diamond is the poorest 
conductor. In the 2D limit, where G oc l/\/c, the differences between the materials are not as 
strong as in the 3D limit. One could say that all materials start to resemble each other in the 
ultrathin range d < 10 nm. 

On the other hand, in the 3D limit, G increases linearly with d, the common sense result. 
Therefore, we could try to define a materials parameter G/ld, the conductance per unit emission 
area, for all thicknesses. This is plotted in Fig. [2] (b). As we see, in the 3D limit G/ld is not 
a function of d, and is truly a materials parameter. However, when the membrane becomes 
thinner, G/ld starts increasing strongly. This means that thin membranes are actually more 
effective thermal conductors, per unit emission area. This increase in ballistic "conductivity" is 
striking: for diamond it can be four orders of magnitude from d = l/xm to d = 1 nm. Therefore, 
a way of increasing the ballistic thermal conductance in a given volume would be to use a 
structure composed of many parallel thin membranes of a material with low speed of sound. 

The dependence of G on temperature is plotted in Fig. [3] for d = 100 nm. In the 3D limit, 
G is naturally proportional to T'^. In the 2D limit however, the power of the temperature 
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Figure 3. Ballistic conductance G as 
function of temperature. In the 3D limit, 
G oc T3, while in the 2D limit we have 
G oc T3/2. 
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Figure 4. Logarithmic derivative of the 
ballistic conductance a = dln{G) / d\n{T) as 
function of temperature. In the 3D limit 
a = 3, while in the 2D limit a 3/2. 



dependence does not just drop by one, but instead we have G oc T^/^, and the crossover 
temperature is between 50 .. 400 mK for d = 100 nm. To give a better picture of how the 
power of the temperature dependence of G changes with temperature, we plot the logarithmic 
derivative dln{G) / dln{T) of G in Fig. d 

5. Conclusions 

For a given temperature, ballistic thermal conduction is maximized by choosing a material with 
low speed of sound, and by using several parallel thin membranes, instead of a bulk slab of the 
same total thickness. A bulk slab of diamond of thickness 30 iim has the same conductance as a 
2 nm thick lead membrane. The lead membrane is thus about 15 000 more effective conductor. 
Applications in cooling ultrasensitive thermal detectors can be envisioned. 
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